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(54) Rapid high throughput spectrometer and method 



(57) A fluorescence spectrometer comprises a laser 
(12) and at least one beam splitter (18) positioned to 
receive a light beam from the laser and to divide it into 
several first light beam portions (16). Dichroic mirrors 
(17) are positioned to separately receive the first light 
beam portions and to reflect the beam portions at an 
angle to the first light beam portions. The first light beam 
portions respectively enter different spectometer units 
(20a-20d) for the irradiation of separate samples and 
sensing fluorescence emitted from the samples. Trans- 
parent chambers are provided for holding the samples. 
Objective lens systems are respectively positioned in 
the path of the reflected beam portions to respectively 
focus each reflected beam portion to a point within one 
of the separate transparent chambers. Lenses are po- 
sitioned to receive fluorescence from a sample for test- 
ing within the transparent chambers and to respectively 
focus the fluorescence at pin holes in opaque partitions. 
The lenses are positioned to receive the fluorescence, 
which passes back through the objective lens system 
and the dichroic mirror. Light detectors are each respec- 
tively positioned adjacent to one of the partitions, with 
one of the partitions respectively positioned between 
each of the lens and the light detectors to permit each 
light detector to sense fluorescence through the pin 
hole. Electronics are provided to receive and process 
signals from each light detector. Structure is provided to 
permit high speed data collection from a large number 
of samples in separate, transparent chambers. 
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Description 

BACKGROUND OF THE INVENTION 

5 [0001] The invention relates to improvements in a multiple channel spectrometer capable of quickly analyzing large 
volumes of samples by the fluorescence emitted bythesamples. While many different systems for causing fluorescence 
emission and processing of the fluorescence may be used, one currently preferred method is known as fluorescence 
correlation spectroscopy (also known as fluorescence fluctuation spectroscopy). Examples of this technique are as 
described in Schrof et al. U.S. Patent No. 5,815,262. Other examples of fluorescence cross-correlation spectroscopy 

10 are shown in the articles by Andre Koltermann et al. from the Proceedings of the National Academy of Science , Volume 
95, pages 1421-1426, February 1998 entitled "Rapid Assay Processing by Integration of Dual Color Fluorescence 
Cross-Correlation Spectroscopy: High Throughput Screening for Enzyme Activity"; and the article by Petra Schwille 
et al. from the Biophysical Journal , Volume 72, pages 1878-1886, April 1997, entitled "Dual Color Fluorescence 
Cross-Correlation spectroscopy for Multi Component Diffusional Analysis Solution". 

15 [0002] In the Schrof et al. Patent No. 5,81 5, 262, fluorescence correlation spectroscopy (FCS) is described in which 
an excitation laser beam passes through a series of wells (sample chambers) arranged in linear array. The beam is 
refocused prior to entering each sample chamber, so that fluorescence takes place in a very small area, which area 
is monitored and sensed to obtain fluorescence data simultaneously from a plurality of samples. 
[0003] With such a technique, difficulties have been found in the refocusing of light after the laser beam has passed 

20 through the first well or sample chamber. Also, possible change of the laser light beam may take place as it passes 
through the various, separate samples, which may effect the data obtained from particularly the "downstream" samples 
in the linear array. Also, two photon excitation, as is used in the above cited patent, generally requires a very expensive 
laser, while one photon excitation permits the use of cheaper lasers. 

[0004] Additionally most of the prior art systems utilize a conventional hardware correlator to receive the raw data 
25 from multiwell plates in fluorescence correlation spectroscopy. Data from the hardware correlator is passed to the 
computer processing unit for determination of the parameters of interest, which may include diffusion coefficients and 
the concentrations of the components. However, several limitations hinder the use of a hardware correlator in this 
manner. Specifically, when a particulate or an aggregate is present in the solution being analyzed, and it passes through 
the tiny region illuminated by the focused laser beam, the entire calculated autocorrelation function resulting from the 
30 data has to be rejected, as it is altered by the temporary presence of the particulate or aggregate. As the result, the 
measurement on that specific well from where the data comes has to be rejected, and, more likely, measurement of 
the entire multiwell plate has to be acquired again. 

[0005] Also, the raw data may contain useful information and features that are completely and definitively lost once 
the data has been processed through the hardware correlator. By system of this invention, the user can examine the 
35 raw data after they have been acquired. Specifically, it may turn out that higher order correlation functions may be of 
more interest in describing the molecular interactions occurring in the sample solutions being analyzed through the 
spectrometer. Thus, if one can keep the raw data, this permits the user to further analyze it with more complicated 
analysis models when that is desired. 

[0006] The analysis can be easily implemented in an automated fashion by properly designed software for a spec- 

40 trometer, particularly for high throughput screening instruments. 

[0007] Also, the first order autocorrelation function of fluorescence correlation spectroscopy is typically determined 
by using the time-mode, which is the traditional way of calculating the function. See particularly Thompson, Fluores- 
cence Correlation Spectroscopy, in "Topics in Fluorescence Spectroscopy" , Volume 1 (J.R. Lakowicz, Editor), Plenum 
Press, New York 1 991 , pages 337-41 0. Time mode operations tend to limitthe precision in determining concentrations, 

45 and increase the data acquisition time. 

[0008] It is desired to perform drug and other screening at very large rates of analysis. For example, current tech- 
niques can allow screening up to 50,000 to 1 00,000 compounds a day. However, it would be desirable in the field of 
high throughput screening to significantly increase the capacity of spectrometers to process large numbers of com- 
pounds. Many of the current high throughput screening apparatus are manufactured by L. J.L. Biosystems of Sunnyvale, 

50 California; Aurora Biosciences of San Diego, California; Molecular Devices of Sunnyvale, California; and Packard 
Instruments of Meriden, Connecticut. 

[0009] Fluorescence emission is usually preferred for such high throughput uses due to the overall sensitivity when 
compared to other techniques such as absorption measurements. Another advantage of using fluorescence is the 
detection method of availability of a range of fluorophores that can be used as extrinsic probes. Typically, five param- 
55 eterscan be measured when using afluorescence technique, namely the intensity of the excitation spectra, the intensity 
of the emission spectra (each at selected wavelengths), the polarization of the excitation spectrum, the quantum yield 
of fluorescence, and the decay time of the excited level. 

[0010] Fluorescence correlation spectroscopy was originally proposed by Magdeet al., Thermodynamic Fluctuations 
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in a Reacting System: Measurements by Fluorescence Correlation Spectroscopy, Physical Review Letters , Volume 29 
(1972), pages 705-708. In this technique, the temporal fluctuations of the detected fluorescence signal (that is time- 
dependent, spontaneous intensity fluctuations of the fluorescence signal in the typically tiny observation volume) are 
detected and analyzed to obtain information about the processes occurring on a molecular scale. These intensity 

5 fluctuations and the volume under observation may arise from Brownian motion, flow, and chemical reactions. During 
the past years, fluorescence correlation spectroscopy has been utilized to measure transitional diffusion coefficients, 
rotational diffusion coefficients, kinetic rate constants, molecular aggregation, and molecular weights. An article by 
Thompson et al., presents a review of the technique (N.L. Thompson et al., Fluorescence Correlation Spectroscopy, 
in " Topics in Fluorescence Spectroscopy ", Volume 1 (J.R. Lakowicz, Editor Pleanum Press, New York 1991, pages 

10 337-410). 

DESCRIPTION OF THE INVENTION 

[0011] By this invention, an apparatus and method are provided for carrying out high throughput screening of active 
15 compounds, typically using fluorescence correlation spectroscopy, although other techniques may be utilized making 
use of this invention. A fluorescence probe is excited through a one photon or multi photon excitation process. The 
light source may be lamp such as an xenon arc or deuterium lamp, or a laser such as continuous wave lasers: i.e., 
argon-ion, krypton-ion, helium-neon, helium-cadmium, or other lasers. Pulsed lasers may also be used such as nitrogen 
lasers or mode-locked lasers, diode lasers, or lasers placed in an array. In each of the possible radiation sources, the 
20 light source should be capable of delivering radiation at a particular wavelength or wavelengths that excite the fluo- 
rescence probe through one photon or multi photon excitation processes. Typically, such excitation wavelengths may 
range from 200 nm. to 5,000 nm. 

[0012] By this invention, a fluorescence spectrometer is provided which comprises a laser; and at least one beam 
splitter (which may be a prism-type beam splitter, a fiber optic system, or similar device for accomplishing the beam 
25 splitting function) positioned to receive a light beam from the laser and to divide the beam into a plurality of separate, 
first light portions. Thus, multiple first light beam portions are provided, typically spaced from and extending parallel to 
each other. 

[0013] Dichroic mirrors are positioned to separately receive the first light beam portions and to reflect the beam 
portions at an angle to the first light beam portions, typically perpendicular. A plurality of transparent wall chambers 

30 (such as wells) are provided for holding samples to be analyzed. Typically, a conventional multiple well sample plate 
may be used. Objective lens systems are provided (typically for confocal use), for example microscopes, which systems 
are respectively positioned in the path of each of the reflected first beam portions, to respectively focus each reflected 
beam portion to a point within one of the separate, transparent chambers, to elicit a fluorescent response from the 
sample for testing in the chamber. 

35 [0014] Lenses are respectively positioned to receive the fluorescence from the sample for testing (which samples 
are within the transparent chambers), and to respectively focus the fluorescence at pinholes in respective, opaque 
partitions. The lenses are positioned to receive the fluorescence which passes back through the objective lens system 
and the dichroic mirror, and then to focus the fluorescence at a respective pinhole. The dichroic mirror is selected so 
that it is transparent to at least one wavelength of the fluorescence. 

40 [0015] Light detectors are each respectively positioned adjacent to one of the pinholes, with that pinhole being re- 
spectively positioned between the lens and the light detector of the particular system, to permit each light detector to 
sense fluorescence through the pinhole. Electronics are then provided to receive and process signals from each light 
detector. 

[0016] The spectrometer of this invention will comprise a plurality of such individual systems, typically four or eight, 
45 for simultaneous processing of sample in separate chambers or wells, to contribute to a high throughput system. 

[0017] Typically, a plate comprises the transparent chambers described above, plus a plurality of other transparent 
chambers for holding samples, typically carried by the same plate. An x-y movement device carries the plate and 
permits a first group of the chambers to be respectively and simultaneously exposed to the first light beam portions for 
analysis of the resulting fluorescence. Then, the x-y device may be moved in a two dimensional, planar manner for 
50 exposure of another group of chambers of the plate to the first light beam portions, for their analysis. 

[0018] Preferably, the system is set up so that fluorescence from each transparent chamber which passes through 
the objective lens system and the dichroic mirror extends in a straight line, passing further through the objective lens 
system and the pinhole to the light detector in the same straight line. 

[0019] Accordingly, an array of such detector systems, comprising the above items to control and process multiple 
55 first beam portions, can be used to systematically and simultaneously analyze the contents of multiple chambers of a 
plate, followed by translation of an x-y table which carries the plate, for analysis of another set of the chambers, until 
the multiple chamber plate (up to 384 or more such chambers) has been completely and rapidly analyzed. By way of 
advantage, each of the first beam portions produced by this invention pass through only one chamber and sample per 
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individual analysis, so there is no possibility of chamber data cross contamination resulting from a light beam that 
passes through multiple chambers to elicit a fluorescent response, as in the patent cited above. 
[0020] Preferably, the light excitation and detection scheme used in this invention follows the confocal design of 
microscopy. Data acquisition can be achieved by measuring the time intervals between the photons reaching the de- 

5 tector and building a histogram of the detected counts. The volume in the sample from which photons are detected 
can be of extremely small volume, for example, about 0.1 to 10 femtoliters, which is on the order of the volume of a 
bacterium. The photon counting histogram provides the concentrations of the molecular species present in the solution 
and the number of photons emitted by each species. The invention can also derive simultaneously the autocorrelation 
function or higher-order correlation functions of the parameters of interest for rapid high throughput. Screening can be 

10 acquired on such small observation volumes over a time of about one to five seconds. Direct measurement of single 
molecules and the kinetics involved over any desired time scale can be achieved by this technique. The lower limit of 
detection is due to impurities and buffer contamination, plus the effect of the particulate resulting from insoluble com- 
pounds. 

[0021] Data acquisition and analysis of the data may comprise the steps of: acquiring a fluorescence measurement 
15 over time from a sample in the form of electronic raw data. This electronic raw data may be picked up by a light detector, 
sent to a preamplifier discriminator, and then sent to a computer, where the raw data may be stored, in accordance 
with this invention, and contrary to the procedure used with the conventional hardware correlator. 
[0022] The electronically stored raw data is then processed using a first algorithm, without erasing the electronically 
stored, original raw data. Then ; as an advantage of and by means of this invention, the stored raw data may be re- 
20 processed using one or more additional algorithms, or the data may be reprocessed with the first algorithm, without 
the need of acquiring a complete set of new data as may be required in the prior art. 

[0023] In fluorescence correlation spectroscopy (FCS), the fluorescence signal F(t) as a function of time is measured 
as the raw data. The temporal autocorrelation of the fluorescence fluctuations, which is a measurement of the average 
temporal duration of the fluorescence fluctuations, is determined. Typically, the normalized autocorrelation function is 
25 defined as: 



[0024] G(D ) decays in time. The rate of the decay and the shape of the curve contain information about the mech- 
35 anisms and the rates of the processes that generate the fluorescence fluctuations. The observed fluctuations of the 

fluorescence signal obey Poisson statistics with the amplitude of the average fluctuation proportional to N 1/2 , where N 

represents the number of molecules in the observation volume. Nanomolar concentrations can be detected. Using 

FCS with two-photon excitation, the observation volume ranges typically from 0.1 to 1.0 x 10 -15 liters (or femtoliter. 

abbreviated as "fl"). This small observation volume allows the direct measurement of single molecules and the kinetics 
40 involved over a time scale extending from hundreds of nanoseconds to seconds or hours. The lower limit of detection 

is due to impurities and buffer contamination and the effect of the particulate generated from insoluble compounds. 

[0025] It is generally preferred to acquirethedata using a photon-modetechnique ratherthan atime modetechnique, 

although either method may be used in accordance with this invention. 

[0026] In the photon mode technique, as is known per se, the detector records the time delay between one photon 
45 and the next photon from the fluorescence arriving to the detector. In this implementation of the data acquisition, the 
"clocks" are the events to be recorded, and the photons are the starts-stops which define each interval. 
[0027] On the other hand, in time mode techniques, the detector counts the number of photons collected from the 
sample in a specific time interval. Typically, the instrument uses about 256 different time intervals. The length of a time 
interval is specified by the user through the software which controls the operation. Thus, the photons are the "events" 
50 to be recorded at each time interval. The "clocks" are the arbitrary starting and ending stops defined by the electronics 
to create each desired time interval. 

[0028] By way of advantage, the method of this invention may be used in conjunction with a one photon excitation 
technique, in which the excited level of the molecule is created by one photon contrary to the situation where two or 
more photons are needed to create the necessary excitation to achieve fluorescence in the sample. One photon ex- 
55 citation techniques can use significantly less expensive lasers. 

[0029] In this invention, a data acquisition card may be provided in the computer used for processing the data, in 
which the data acquisition card acquires the raw data from the light detector and arranges for its storage in the computer 
memory. Then, the computer may calculate the desired autocorrelation function or higher-order correlation functions. 
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This provides the advantage that bad data may be eliminated without losing the entire data set pertinent to the sample 
well under examination, contrary to the situation of a hardware correlator. The bad data may result from the presence 
of a solid particle in the tiny volume being examined for fluorescence. Additionally, the raw data may be rerun after 
acquisition for recalculation based on higher order correlation functions, if that is needed. Such a data acquisition card 
5 can also have the capability of acquiring the data in either the time mode or in the photon mode, as may be desired. 
[0030] Also, as previously discussed, rapid, multiple, simultaneous analysis of sample chambers may take place in 
accordance with this invention for high volume screening operations. 

DESCRIPTION OF DRAWINGS 

10 

[0031] 

In the drawings, Fig. 1 is a schematic view of a fluorescent spectrometer for FCS, shown in a plan view; 
Fig. 2 is a partially schematic bottom plan view showing the four spectrometer units similar to the units of Fig. 1 
15 but in a rectangular, not a linear arrangement, and their relationship with a microwell plate for holding the samples; 

Fig. 3 is an elevational schematic view of a single spectrometer unit, for simultaneous use along with other similar 
spectrometer units in the device of Figs. 1 or 2; and 

Fig. 4 is a flow chart for software processing of data obtained with the devices of Figs. 1 -3. 

20 DESCRIPTION OF SPECIFIC EMBODIMENTS 

[0032] Referring to Fig. 1 , a multi-channel high throughput screening (HTS) spectrometer 1 0 is shown in schematic 
form, used for FCS. 

[0033] A laser 12 serves as a source of excitement radiation for fluorescent samples. Light beam 14 from laser 12 
25 is split into four separate first light beam portions 1 6 by a group of beam splitters 1 8 of conventional design and mirrors 
1 7. The first light beam portions 1 6 respectively enter different spectrometer units 20a, 20b, 20c, 20d for the irradiation 
of separate samples and the sensing of fluorescence emitted from the samples. The detailed structure of the spec- 
trometer units 20a-20d may be identical, one being specifically shown in Fig. 3. 

[0034] Alternatively, the system of beam splitters 18 and angle mirrors 17 may be replaced by a multiple light cable 
30 system, in which light is transmitted from the laser by the cable system to provide individual first light beam portions 
directed to the respective spectrometer units 20a-20d. 

[0035] There is no intrinsic limit on the number of individual spectrometer units 20 which may be irradiated by a single 
laser 12. For example, at least eight separate spectrometer units may be serviced by laser 12 and a light distribution 
system 17, 18 (or light cables), up to about 150 spectrometer units and more if desired, depending upon the power of 
35 the laser. 

[0036] Fig. 2 shows a plan view of a system which is basically identical with Fig. 1 except that the respective, cor- 
responding spectrometer units 20a', 20b', 20c', 20d' are arranged in a square rather then a straight line pattern by 
appropriate modification of the light distribution system 17, 18 (or light cables). The spectrometer units point at the 
transparent bottoms of chambers 22, which are present in a checkerboard array as part of a commercially available 

40 384 well plate 24, so that 384 separate samples may be analyzed by the fluorescence spectrometer of this invention. 
Plate 24 is carried in an open bottom X-Y table 26, which is capable of precise movement in both horizontal directions, 
so that after a typical 1-5 second irradiation of the four wells 22 that are being analyzed by the respective spectrometer 
units 20a'-20d', X-Y table 26 can be horizontally shifted so that four new wells 22 may then be tested by the four 
spectrometer units. Accordingly, the system of this invention is capable of high volume production of spectrograph ic 

45 analysis of a large number of samples. X-Y tables are well-known, commercially available devices. 

[0037] Referring to Fig. 3, an elevational view of the various spectrometer units 20, and their design, may be seen. 
Laser 12 emits light beam 14, which light beam is repeatedly split by a light distribution system comprising a series of 
beam splitters 18, to form the respective separate, first light beam portions 16, as in Fig. 1 . By the appropriate arrange- 
ment of beam splitters 18 [and diagonal mirrors 17, as desired], first light beam portions 16 may be directed in any 

50 pattern desired, linear, rectangular, or otherwise for transfer into spectrometer units 20. 

[0038] As shown, one of the beam portions 1 6 enters unit 20 and strikes dichroic mirror 28, being reflected upwardly 
as a beam portion 1 6a. Dichroic mirror 28 is a commercially available mirror having the characteristic of being reflective 
to light at some wavelengths and transmissive to light at other wavelengths. Of course, first light beam portion 16 is 
selected to be of a wavelength that is reflected by dichroic mirror 28. 

55 [0039] Light beam 16a passes through objective 30, which generally functions as a microscope, focusing the light 
beam to a focal point 32 residing within microwell 22 and within the sample contained therein. As previously stated, 
the volume at focal point 32 can cover a very small volume of material, on the order of a one femtoliter volume approx- 
imating the volume of a bacterium. The laser intensity is adjusted so that photo bleaching of the sample is avoided, 
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resulting from excessive beam intensity. The optimum beam intensity depends upon the sample concentration plus 
the sample quantum yield, which is the number of photons emitted by the sample divided by the number of photons 
passed into the sample. Some fluorophores may have a very low quantum yield, which naturally will then require a 
larger laser beam intensity. 

5 [0040] The fluorescent response may take place from one or more of the molecules within the tiny, focused volume 
32, the fluorescence being at a characteristic wavelength that is typically different from the wavelength of the laser 
beam 16,1 6a. Some of this fluorescence passes back through objective lens 30 and through dichroic mirror 28 without 
much reflection, since it is of a different wavelength. The dichroic mirror is of a type selected to be transparent to the 
particular fluorescent wavelength. It is desirable for the fluorescence to pass back through the same objective lens 30 

10 that the irradiating radiation passes through in the other direction, since when a second objective lens is used, it would 
be difficult for the second objective lens to align and focus upon the same volume 32 that the first objective lens 30 is 
focusing on. 

[0041] After passing through the dichroic mirror 28, filter 34 blocks light of wavelengths other than the desired wave- 
lengths. The fluorescent light is focused by lens 36, so that it passes through a pinhole 38 in an opaque partition 40 

15 at or very near the focal point of lens 36 under the particular conditions. 

[0042] The amount of the volume of sample probed by the instrument (typically at the femtoliter range) is determined 
by the optics of the system, plus the diameter of pin hole 38. If desired, an instrument incorporating this invention may 
carry a series of different pinholes which may be alternatively placed into the position shown in Fig. 3, to provide a 
variation of the volume of sample 22 thus probed by spectrometer unit 20. 

20 [0043] A light detector 42 of conventional design is monitoring pinhole 38 from the other side of partition 40. When 
light is detected, it will be typically only fluorescence from the sample, since other wavelengths of light are removed 
by filter 34. Signals from the light detector 42 pass through cable or wire 44 to preamplifier discriminator 46, and the 
signals from discriminator 46 pass through cable or wire 48 to computer 50. In the computer, the fluorescent signal, 
measured as a function of time, is processed to obtain a temporal autocorrelation of the fluorescence fluctuations, 

25 which is a measurement of the average temporal duration of the fluorescence fluctuations. Using the well-known nor- 
malized autocorrelation function, it becomes possible to directly measure fluorescence characteristics of single mole- 
cules, and the kinetics involved over a time scale extending from hundreds of nanoseconds to seconds or hours. 
[0044] Referring to Fig. 4, the data is processed by means of the above technique as follows: the measurement 60 
is acquired on a well 22 in the manner previously described. The data passes as a data stream into computer 50 as 

30 at 62. A math algorithm 64 examines the raw data for the presence of large signal variations indicating contaminants. 
If present, large variations in the signal are removed by a software filter 66, and the data passes through math algorithm 
64 again. If no large signal variations are found, a math algorithm 68 calculates the autocorrelation function of the 
acquired data. Parameters of interest are determined, such as diffusion, coefficients, concentrations, binding constants, 
etc. through software 70. 

35 [0045] At decision point 72 of the flowchart, results are either accepted in accordance with predetermined criteria of 
algorithm 73. If the criteria are not achieved, the math algorithm 74 calculates higher order correlation functions through 
the computer software. Again a decision point 76 is reached. If the results are now in compliance with the criteria of 
the algorithm 78, the results are accepted and provided to the user. Otherwise, the algorithm may be predetermined 
to provide another analysis model 80 for the acquired data. This recursive loop of decisions including decision point 

40 72, 76, 82 may continue for as long as may be desired, until the results fit the desired criteria or there are no other 
algorithms or analysis techniques available in the program to process the data. 

[0046] By such techniques, a wide variety of studies of the behavior of matter, living and otherwise, can be performed. 
The translational mobility of cytoplasm in living cells on latex beads can be determined. Additionally, the presence of 
chemical complexes and their proportional concentration, compared with the separate reactants can be measured. 
45 Drugs and other chemicals can be screened for their presence or absence. Also, a wide variety of theoretical studies 
for physical chemistry and the like can be performed. 

[0047] Specifically, the use of the above apparatus and method may be optimized by the algorithms discussed below: 
[0048] In high throughput screening (HTS) applications important parameters are the sensitivity and the capability 
to separate productive species from background signal and unproductive species fast and efficiently. Until recently. 

50 most of the methods employed either were slow or have very low efficiency. The analysis of fluorescence from single 
molecules could open a new avenue for sensitivity and selectively. The major concept in this field is how to analyze 
signals from different molecular entities. If the molecules of interest are fluorescent, the methods of single molecule 
detection and analysis can provide a powerful new way to perform an analysis of the sample on a molecule-per- 
molecule basis. To achieve this goal, 1) we must restrict the analysis volume to the conditions that one or at most few 

55 molecules be present in that volume, 2) we must collect very efficiently the fluorescence emitted from that small volume, 
3) we must devise fast on-line algorithms to analyze the fluorescence signals and sort the results according to the 
different molecular populations present in the sample. Since the fluorescence properties are analyzed on the flight, 
there is in principle no limit on the number of different fluorescence population that we can separate: the actual limit is 
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given by the signal-to-noise ratio. For example, using this principle of molecular sorting, we can separate one bright 
molecule from a large number of dim fluorescent molecules even if the concentration of the bright molecule is very low 
compared to the dim molecules. It is precisely this principle that we are applying to the HTS apparatus to increase the 
sensitivity and the speed of data analysis. In reality, it is not always possible to reduce the concentration of the fluo- 
5 rescent molecules to the limit in which only one molecule is present in the analysis volume and the analysis of the 
statistics of the observed fluorescence signal requires complex mathematical methods. However, we were able to 
reduce the mathematical complexity to simple formulas that can be handled by modem computers to allow on-line 
analysis. In the following we describe how we achieve the 3 conditions outlined above for this new method to work 

10 Small volume. 

[0049] We have exploited the confocal principle and/or multiphoton excitation to achieve fluorescence excitation 
volume which are diffraction limited. The volumes that can be obtained using current state-of-the-art microscopy ob- 
jectives and common lasers are on the order of 0.1fl_. A volume of that size contains on the average 0.06 molecules 
15 for a sample at a concentration of 1 nM. These concentrations are typical in HTS applications. 

Efficient light collection. 

[0050] We are using relatively high numerical aperture objectives and a particular cell for HTS that allow us to use 
20 short working distance objective. Also the light detector we use (avalanche photo diode) has very high quantum effi- 
ciency in the spectral region of most dyes used in HTS. We also use relatively high excitation power to saturate the 
absorption of the dye. Absorption saturation is easily obtained using conventional lasers both in one or multi-photon 
excitation. As a consequence we can routinely obtain counting rates exceeding 100000 counts/sec/molecule. During 
the fluorescent burst when a dye molecule is under the laser beam which last on the order of 1 ms or longer, assuming 
25 that the dye is attached to a macromolecule, we collect on the average about 50-1 00 photons. This number of photons 
collected is sufficient to provide statistics for the separation of dye molecule which differ in brightness by about a factor 
of 2. 



30 



Fast algorithms 

[0051] We have developed a series of fast mathematical algorithms based on the principle of the analysis of the 
distributions of photon count histograms (PCH) which can analyze on line complex situations in which more than one 
molecule present in the excitation volume. 

35 Additional considerations. 

[0052] The apparatus we assembled allows the simultaneous measurement of the fluorescence signals using two 
light detectors. Therefore, it is possible in principle to measure two separate spectroscopic properties such as differ- 
ential wavelength emission and polarization of a molecule. The principles described below can very well adapt to a 
40 multiple detector system 

Description of the theory 

[0053] Fluctuation experiments observe a physical process that is stochastic in nature. The time series of such a 
45 random signal can be rather complex. Consequently, an analysis based on statistical methods is required to examine 
stochastic processes (Gardiner, 1 985). Traditionally, FCS experiments determine the autocorrelation function g{%) from 
the sequence of photon counts in order to examine the statistics of the time-dependent decay of the fluorescence 
intensity fluctuations to their equilibrium value. The determination of the autocorrelation function from the raw data is 
a data reduction technique. Hence, some information encoded in the time series of the photon counts is lost. While 
50 the autocorrelation approach is the method of choice for characterizing kinetic processes embedded in the stochastic 
signal, it lacks information regarding the amplitude distribution of the intensities (Bendat and Piersol, 1971). Here we 
consider a different data analysis approach based on the amplitude distribution of the fluorescent intensities. Experi- 
mentally photon counts rather than intensities are measured and the statistics of the photon count amplitudes must be 
considered. We are specifically interested in the probability to observe k photon counts per sampling time. This prob- 
55 ability distribution of photon counts is experimentally determined by the photon counting histogram (PCH). The analysis 
of fluorescence fluctuation experiments by PCH had been recently introduced by our group (Chen et al., 1 999). In this 
contribution, we focus on the theoretical foundation of the photon counting histogram and its experimental realization. 
We develop the theory for a single fluorescent species and confirm it experimentally. In addition we apply the PCH 
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method to separate heterogeneous sample compositions. Therefore, the theory is generalized to include multiple spe- 
cies. The autocorrelation function g(%) describes the fluctuations in the time domain, but with the exception of g(0) 
lacks amplitude information. The photon counting histogram, on the other hand, characterizes the amplitude distribution 
of thefluctuations, but lacks kinetic information. Thus, autocorrelation and PCH analysis are complementary techniques 
and PCH should be able to separate a mixture of species based on a brightness difference between the components 
regardless of their diffusion coefficients. We discuss the use of PCH to separate a mixture of components and consider 
sample conditions, such as the molecular brightness and the particle concentration. We demonstrate the technique by 
resolving a binary dye mixture from the photon count distribution. Furthermore, we apply PCH analysis to a mixture of 
biomolecules with either one or two fluorescent labels attached and resolve them experimentally. 

Theory 



[0054] To describe the PCH of a freely diffusing species we first consider a single, diffusing particle enclosed in a 
small box of volume V 0 . The PCH of a single particle depends explicitly on the beam profile of the excitation light. In 
15 the next step, several particles are added to the box and the corresponding PCH function determined. To describe an 
open system with Poissonian number fluctuations, the boundary condition of the volume is removed and particles are 
allowed to enter and leave the box. We expanded the model to cover multiple species and calculate the PCH for a 
number of beam profiles. However, first we consider the statistics of the photon detection process. 

20 PCH and the theory of photon detection 

[0055] Fluorescence fluctuation experiments are typically performed using single photon counting techniques. Each 
detected photon gives rise to an electric pulse, which is sent to the data acquisition system. Thus, the fluorescence 
intensity reaching the photodetector is converted into photon counts. The photon detection process changes the sta- 
25 tistics of the measured intensities by adding shot noise to the photon counts (Saleh, 1978). For example, a constant 
light intensity at the detector l D gives rise to a Poissonian distribution of photon counts k, 

PoK*.<*»=^ • (D 

[0056] Shot noise is a random Poisson point process, and reflects the discreteness and statistical independence of 
the photoelectric detection process (Snyder 1975). The factor i\ f is the proportionality factor between the average 
photon counts {k) and the constant intensity l D at the detector, (k) = j\ ( I d . It incorporates the detection efficiency and 
35 the sampling time interval. The variance 



40 

of a distribution serves as an indicator of its width. Forthe Poisson distribution, the mean and the variance are equal, 



[0057] A Poisson distribution describes the photon count statistics for light of constant intensity. Fluctuations in the 
light intensity change the photon counting statistics and the corresponding distribution of photon counts was first de- 
scribed by Mandel (Mandel, 1958), 



CO 

p(*) = JPoi(i 5 tj, I D ) p(I D ) dI D . (2) 

0 

55 

[0058] The probability of observing k photoelectron events is given by the Poisson transformation of the intensity 
probability distribution p(/ D ). However, to be more precise, photon detection involves a short, but finite sampling time 
At s , which along with the detector area A needs to be integrated over to yield the energy W, 
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t+&t s 

W{t)= J \l(r,t)dkdt. (3) 

( A 

5 

[0059] Thus, technically we observe energy fluctuations of the collected light, rather than intensity fluctuations. How- 
ever, if the time scale of the intensity fluctuations is longer than the sampling time At s , then the energy fluctuations 
track the intensity variations. We will assume for the rest of the paper that the sampling time At s is chosen fast enough, 
so that the energy fluctuations track the intensity fluctuations of interest. We also assume a stationary process, so that 
10 there is no explicit time dependence to the statistical properties of the photon counts, and a detector area A small 
enough, so that the intensity field is essentially constant across the detector surface. 

[0060] Mandel's formula (eq. 2) essentially describes a superposition of Poissonian distributions scaled by the re- 
spective probability of the intensity distribution function p(/ D ). The photon count distribution p(k) is now characterized 
by a variance 

15 

K) 

greater than its mean value, 

20 

(Ak 2 )>(k), 

which is classified as super-Poissonian (Teich and Saleh, 1988). Thus, intensity fluctuations lead to a broadening of 
25 the photon count distribution with respect to a pure Poisson distribution (fig. 1). The variance 

30 of the intensity distribution p(/ D ) determines the variance of the photon counts 

35 (Mehta, 1970). Thus, as the strength of the intensity fluctuations increases, so does the broadening of the photon 
counting histogram. The changes in the shape of the histogram from a Poisson distribution are characteristic for the 
intensity distribution p(/ D ) and it is possible to infer some characteristic properties of the light source from the photon 
counting histogram. This approach has been used in the past to investigate the scattering of light (Bertolotti, 1973), 
the twinkling of stars (Jakeman et al., 1978) and the fluctuations of laser light (Risken, 1970). Here we follow the same 

40 approach and directly model the photon counting histogram to describe fluorescence fluctuation experiments of freely 
diffusing particles. 

[0061] The experimental setup typically involves a microscope and an objective, that is used to focus the excitation 
light in order to achieve a small spatial volume. A detector collects the fluorescent light emerging from the excitation 
volume. The effect of the microscope optics and the detector on the shape of the observation volume is characterized 
45 by the point-spread function (PSF) of the instrument. The shape of the PSF influences the photon count distribution 
and will be considered for a few PSF conventionally used in confocal and two-photon detection (Qian and Elson, 1 991 ; 
Berland et al., 1995). In our context, it is more convenient to define a scaled PSF, PSF, such that the volume of the 
scaled PSF, V PSF = JPSF(r)dr, is equal to the volume defined for FCS experiments (Thompson, 1991). 

50 PCH of a single particle 

[0062] We consider a single particle diffusing in an enclosed box with a volume V 0 large enough, so that it essentially 
contains the PSF (fig. 2a). Based on Mandel's formula we derived an equation which expresses the photon count 
distribution as a function of the PSF (Chen et al., 1 999), 

55 
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p 0) (k;V 0 ,s) = |Poi(* 5 rf>SF(r))p(r)dF = i- |Poi(*, sT>SF(r))dr , (4) 

5 

where p(r) describes the probability to find the particle at position ?. For a freely diffusing particle, the probability equals 
1/ V 0 inside the box and zero outside of the box. The photon counting probability p( 1 )(/c; Vq,z) for a single particle enclosed 
in a volume V 0 depends on the shape of the PSF and the parameter e. The meaning of this parameter is best illustrated 
by considering the average photon counts (k) of the PCH p( 1 )(/c; t/ 0 ,e), 

10 

(k) = JL$PSF(r)dr = z*f^. (5) 

15 

[0063] The average photon counts are thus determined by the product of e and the probability to find the molecule 
within the volume of the point spread function V PSF . Therefore e describes the molecular brightness, which determines 
the average number of photon counts received during the sampling time Af s for a particle within the observation volume 
V PSF . The average photon counts received (k) scale linearly with the sampling time. Therefore, the ratio e sec =e/At s is 
20 independent from the somewhat arbitrary sampling time At s . The parameter e sec expresses the molecular brightness 
in photon counts per second per molecule (cpsm) and allows a more convenient comparison between different exper- 
iments. 

PCH of multiple particles 

25 

[0064] Now let us consider N independent and identical particles diffusing inside a box of volume V 0 (fig. 2b). If one 
could follow a particular particle individually, the PCH of this particle would be given by pC)(/r; V 0 ,s) according to eq. 5. 
For N independent particles the corresponding PCH, p( N )(/c; l/ 0 ,e), is given by consecutive convolutions of the single 
particle PCH functions p( 1 )(/c V Qt s) (Feller, 1957), 

30 

y™(k;V 0 ,z) = ( p (l) <8> ••■ ® y m ) (k;V Q ,z) . (6) 

N-times 

35 

PCH of particles with number fluctuations 

[0065] The assumption of a closed system, in which particles diffuse inside a box, does not describe the experimental 
situation, unless the reference volume includes the whole sample. But a macroscopic reference volume would require 
40 the evaluation of an astronomical number of convolutions according to eq. 6. Instead, we choose to consider an open 
system in which particles are allowed to enter and leave a small sub-volume (fig. 2c). The sub-volume is in contact 
with a much larger reservoir volume and the distribution of the number of particles N inside the sub-volume is given 
by a Poisson distribution (Chandrasekhar, 1943), 

p # (A/) = Poi(A/, JV), (7) 

where N describes the average number of molecules within the reference volume V Q . Of course if there are no particles 
in the reference volume, no photon counts are generated and we define the corresponding PCH as, 

50 

pV\k;V 0 ,e-) = Z(k) , with S(*) = |J'*"° . (8) 

55 

[0066] Now we can express the PCH for an open system p(/c; V 0 , N, e) as the expectation value of the N -particle 
PCH p (N) (/e; V 0 ;e) considering Poissonian number statistics, 
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U(k;N PSf ,s) S p(k;V Q ,N ,s) = (p {N) (k;V 0 ,e)} n . (9) 



A. 



[0067] The PCH function p(k; V 0 ,N.e) describes the probability to observe k photon counts per sampling time for an 
open system with an average of N particles inside the reference volume V 0 . The particular choice of the reference 
volume for an open system is irrelevant. It is intuitively clear, that the properties of an open system have to be inde- 
pendent of the arbitrary reference volume V 0 (Chen et al., 1 999). Thus, the photon count distribution should be either 
referenced to an intensive quantity, like the particle concentration, or to some standard volume. We choose the con- 
vention used in FCS, where the volume of the PSF, V PSF , connects the g(0) value of the autocorrelation function to 
the average number of molecules N PSF (Thompson, 1 991). Consequently we drop the V 0 parameter dependence for 
the PCH of an open system and declare a new function U(k; N PSF , e), which characterizes the PCH of an open system 
referenced to the volume of the PSF. The average number of photon counts (k) can be calculated from eq. 9 and is 
given by the product of the molecular brightness s and the average number of particles A/ PSF inside the PSF volume, 



(k) = eN PSF . (10) 



PCH of multiple species 



[0068] So far, only identical particles have been treated. Often more than one type of particle is present in the sample. 
It is straightforward to expand the theory under the assumption that the particles are non-interacting. Let us consider 
the case of two different species for simplicity. If we could distinguish the photon counts emerging from each species, 
we could directly determine the PCH of each species, U(k, , e.,) and n/c; N 2 > ^2- However, we cannot distinguish 
the origin of the photon counts. But as long as the photon emission of both species is statistically independent, the 
PCH of the mixture is given by the convolution of the photon count distributions of species 1 with the one of species 2, 

n^A^/v^e^) = n(/c;A/ 1 ,8 1 )®n(/c;A/ 2 ,e 2 ). (11) 

[0069] For more than two species all single species photon counting distributions are convoluted successively to 
yield the photon count distribution of the mixture. 

PCH for different PSFs 

[0070] The photon count distribution depends on the PSF. Here we report the PCH of a single particle 



P m {k;V 0 ,e) 



for a few PSF of interest. The Gaussian-Lorentzian squared PSF has been used to describe the two-photon excitation 
beam profile for our experimental conditions (Berland et al., 1995), 



PSF 2 gl(a*)= 2 f? x ex P 

7t CO \Z) 



(12) 



[0071] The PSF is expressed in cylindrical coordinates and the excitation profile has a beam waist co 0 . The inverse 
of the Lorentzian along the optical axis for an excitation wavelength of X is given by, 



<s) 2 {z)=a> 2 c l + l— ,withz*=^-. (13) 




[0072] To calculate the PCH of a single particle for a reference volume V 0 , eq. 12 is inserted into eq. 4, but integrated 
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over all space. Integrating over all space is mathematically convenient and ensures the correct PCH for the open 
volume case, since from a mathematical point of view the PSF extends to infinity. However, the PCH of a closed volume 
is only approximately determined. The quality of the approximation depends on the size of the reference volume V 0 . 
If the volume is chosen so, that the contribution of the PSF to the photon counts is negligible outside of the reference 
volume, then the deviation between both functions is small. However, the PCH of a closed volume is from a practical 
point of view only of minor interest, since the PCH of an open volume describes the experimental situation. We refer 
the interested reader to a more detailed discussion of this point by Chen et al. (Chen et al., 1999). The PCH of a single 
particle is then determined for k> 0 by a one-dimensional integral 



12 4 oc 
7t CO. 



V 0 2X 



4s 



\dx , for k>0. 



(14) 



[0073] The integral, which contains the incomplete gamma function y (Abramowitz and Stegun, 1 965), can be eval- 
uated numerically. 

[0074] A second important PSF is the three dimensional Gaussian PSF, which is used extensively to describe con- 
focal detection (Qian and Elson, 1991; Rigleret al., 1993a), 



PSF3dg (x, y, z) = J (*'- y '*) = exp 

-'o 



2(x 2 +y 2 ) 2z 7 



(OX 



(15) 



with an effective beam waist z 0 in the axial direction. The PCH of a single particle is determined in the same way as 
for the other PSF and we derive for k > 0 an expression in the form of a one-dimensional integral, 



pfU^V 0 ,s) = jr^^) r(k,se-^)dx ,for *>0. 



(16) 



[0075] For the measurement of surface processes the PS F is typically approximated by a two-dimensional Gaussian , 



PSF2DG (x, y) = = exp 



2{x 2 +y 2 ) 



(17) 



and the corresponding PCH is given by, 



1 7ZCO 



p«j (k;A 0 ,s) = -^-f- r (k, s) , for k > 0 . 



A 0 2k\ 



(18) 



[0076] Here we reference of course to an area A 0 instead of to a volume. Finally, if the PSF is uniform, then the PCH 
is simply a Poisson distribution, P^ (k;V PSF ,e) = Poi(/c,£). 

Describing PCH with the moment generating function 

[0077] The PCH of a single particle has just been described for photon counts k> 0. To evaluate the PCH for k = 0 
one must determine the sum over all photon counts k> 0 and subtract it from the area of the probability distribution, 
which is normalized to one, 



p (,) (o)=i-Z;. 1 p< i) W- 

[0078] Furthermore, the computation of the PCH for each photon count requires a numerical integration. 
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[0079] We now introduce an alternative formalism to determine the probability distribution of photon counts using a 
moment generating function. This approach not only determines the photon count probability for zero photon counts 
directly, but also expresses the single particle PCH in form of an analytical expression. The description is based on a 
formal relationship between the Laplace and Poisson transform. The Laplace transform L and the Poisson transform 
P are formally defined as, 



P(«)=P|/(jc)]={^/(x)^1, 

J n\ (19) 

F(*)=L[/(x)]= \dxf{xy** . 



[0080] The Laplace transform F(s) of the probability function 1{x) is also its moment generating function (Saleh, 1978), 



(/">-(-!) 



ds" 



(20) 



5=0 



[0081] The Poisson transform can be expressed in terms of the Laplace transform F(s) through the following rela- 
tionship (Saleh, 1978), 



n\ ds" 



(21) 



[0082] Often the Poisson transform can be calculated for n = 0 or 1 in closed form. For example, for the squared 
Gaussian-Lorentzian PSF. we get the following result for n = 1 , 



P2GL(l)=^2F2 



1/4, 3/4 
1/2, 2 

with 



jSS 

K 



V 



8F{s) 



5=1 



ds 



(22) 



[0083] The properties of the generalized hyp ergeo metric function p F q are described in the literature (Slater, 1966; 
Luke, 1 969). The zeroth moment of a probability distribution is equal to one, therefore the moment generating function 
satisfies the condition . F(0) = 1 , according to eq. 20. We can determinethe moment generating function F(s)from eq. 22, 



F(5) = l + r| 



2 F 2 



-3/4,-1/4 
•1/2, 1 



TV 



(23) 



[0084] The photon counting histogram can now be determined from the moment generating function according to 
eq. 21 , using the analytical form of the derivatives of the generalized hypergeometric function and the notation (a) k 
for the Pochhammer function, 
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p2l(*)- 



l + r- 



1- 2 F 2 



'-3/4,-1/4 
-1/2, 1 



,2k 



71 



s k * (-3/4), (-1/4), 



^" r 3(-l/2) t (l), 



,forA: = 0 



2*2 



'-3/4 + *,-l/4 + * 
-1/2 + *, 1 + * 



7t 



(24) 



, for A: > 0 . 



[0085] By using the moment generating function to describe the PCH of a single particle we arrive at an analytical 
solution for p (1) including k= 0 . 

[0086] A second example for the use of the moment generating function is illustrated for the homogeneous PSR In 
this case, a Poisson distribution determines the PCH of a single particle. The PCH for an open system is described by 
a compound Poisson distribution, 



n(*; N,e) = J Poi(£, s N) Poi(iV, N) . 



(25) 



N=0 



[0087] Instead of evaluating the sum in eq. 25 one can use a moment generating function to determine the PCH. 
Here we use the factorial moment generating function 



of the compound Poisson distribution, 



(26) 



[0088] The moment generating function Q(s) allows us to calculate the PCH of the open system analytically by using 
the following relationship, 



(27) 



5=1 



Twofold PCH statistics 

[0089] The photon count distribution characterizes the amplitude fluctuations of the detected photons, but lacks 
kinetic information. We want to demonstrate how to expand the theory of PCH to include time dependence and show 
how PCH is related to the autocorrelation function of FCS experiments. The photon count distribution considered so 
far describes the number of detected photons in a single time interval. However, by studying the joint statistics of photon 
counts from two short time intervals separated by a time delay t, we incorporate time-dependence to the photon count 
statistics (Saleh, 1978). Here we specifically consider the two-fold distribution of photon counts for a diffusing particle. 
The conditional probability that a particle at position r 0 at time f 0 will be found at a later time ^ at position r A is given 
by the solution of the diffusion equation (see for example (Chandrasekhar, 1 943)), 



Pdfa('oMi))=(4*^) 



-3 

2 e 



4Dt 



(28) 



where D is the diffusion coefficient of the particle and x the time difference t A -t 0 . In the following we only consider time 
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lags x> 0 to avoid the complication of shot noise correlations at i = 0. The photon count probability p(k, r 0 ) for a particle 
with molecular brightness e at position r 0 is given by a Poisson distribution, p(kjo) = Poi(k,e PSF(/o)). The average 
photon counts (k(ro)) at position r 0 is proportional to the fluorescence intensity l(r 0 ), (k(r 0 )) = t| /(/q). The propor- 
tionality constant rj describes the detection efficiency and is set to 1 for simplicity. The probability p ( \k 0 ,k^;e,x) to 
observe k 0 photon counts at time f 0 and fc, photon counts at a later time ^ for a diffusing particle with brightness e is 
then determined by (Qian, 1990), 

p (l) (£o,*i;^)= JJdF 0 d^p(rJp^ , (29) 

where p(/o) describes the probability to find the particle at position r 0 . The correlation function of the photon counts 
for a diffusing particle is equal to the correlation function of the fluorescence intensity, 

<* 0 (/)*, (t + r)> = £ £ k 0 k lP W (k 0 , k x ; e, t) = 

jt o =0/t,=0 

JJdW p(Po)/(-o)p d (-o(^^(0)^)= (30) 
(7(/)/(r + r)> 

[0090] The above result was derived for a single particle, but it is straightforward to generalize the result to any 
number of particles or to an open volume with number fluctuations. 

Data analysis 

[0091] We developed an algorithm to calculate the PCH Tl(k) based on the theory presented above and implemented 
it on a computer to fit experimental data. Photon counts are recorded with a home built data acquisition card, which is 
interfaced to a computer. The computer calculates the histogram of the experimental data. The normalized histogram 
represents the experimental photon counting probability distribution p(*)< The statistical uncertainty associated with 
each element of pW is determined by the standard deviation 

The number of data points M collected is typically of the order of 10 6 . The experimental data are fit by minimizing the 
reduced % 2 -function, 

a k (31) 



z 1 - 



k -k -d 

max mm 



[0092] The experimental photon counts range from a minimum value /c mjn , which is typically 0 for most experiments, 
to a maximum number k max , and the number of fitting parameters is given by d. The normalized residuals of the fit are 
determined by 

r(k)=M m=im. 

Single species PCH 

[0093] We tested the theory of the photon counting statistics by comparing it to experiments. The numerical PCH 
algorithm was used to fit the experimental photon counting histograms to the theory as outlined in the data analysis 
section. The experiments were carried out using a two-photon setup and the details of the experiment are described 
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in Chen (Chen, 1999). The measured photon counting distributions agree with the theoretical PCH functions calculated 
forthe experimental setup within the statistical error. The molecular brightness e and the average number of molecules 
shape the photon count distribution in characteristic ways. We first study the influence of the average number of mol- 
ecules N upon the properties of the histogram. 

Influence of the particle concentration 

[0094] The photon count distribution of a fluorescein solution was measured at three different concentrations (fig. 
3). We performed a global fit of all three histograms with e linked, while the average number of particles was allowed 
to vary. The data and the fitted histograms in fig. 3 are in good agreement. The residuals between data and fit for each 
histogram are displayed in units of standard deviation o. The residuals vary randomly and yield a reduced % 2 close to 
1, indicating a good description of the data by the theoretical model. The recovered number of molecules N scales 
with the concentration of the sample. The dashed lines in fig. 3 represent the corresponding Poisson distributions with 
a mean equal to the average photon counts (k) of the sample. The Poisson distribution approximates the PCH for the 
high fluorescein concentration case (fig. 3a). However, lowering the dye concentration to 55 nM (fig. 3b) already shows 
a broadening of the experimental PCH compared to the Poisson distribution, which is clearly visible in the tail of the 
distributions. The deviation of the PCH from a Poisson distribution becomes even more apparent by reducing the 
fluorescein concentration to 5.5 nM (fig. 3c). 

[0095] Thus, the photon counting histogram approaches a Poisson distribution with increasing fluorophore concen- 
tration. This behavior can be readily understood by considering the influence of the molecule concentration on the 
intensity fluctuations. The relative strength of the number fluctuations is given by the ratio between the standard devi- 
ation a and the mean jx of the particle number distribution, 



[0096] The number of molecules inside a small, open volume is Poisson distributed, and the relative strength of the 
particle fluctuations decreases with the inverse square root of the average number of particles N. Thus with increasing 
particle concentration the number distribution approaches a delta function 5(A/- A/). Consequently, the fluctuations in 
intensity become negligible, when compared to the average intensity. The second contribution to the intensity fluctu- 
ations, which is due to the diffusion in an inhomogeneous excitation profile, also vanishes at high particle concentra- 
tions. Any vacancy created by a molecule leaving a position is practically always filled by another molecule moving to 
that position, so that no net change in the fluorescence intensity occurs. Thus, a constant fluorescence intensity dictates 
a Poissonian photon count distribution. 

Influence of molecular brightness 

[0097] We used three different fluorophores, each with its own brightness parameter 8, to illustrate the influence of 
the molecular brightness 8 on the photon count distribution. Each fluorophore sample was made up to approximately 
the same concentration to facilitate the comparison of the different histograms. The count distributions are analyzed 
with the PCH algorithm and are shown together with the fits in fig. 4. Poisson distributions with the same mean as the 
average photon counts are displayed as dashed lines for each histogram. The deviation between the tail of the PCH 
and the Poisson distribution increases with increasing 8. 

[0098] To maximize the deviation between the photon count distribution and the corresponding Poisson function, 
one can either reduce the number of molecules within the excitation volume or increase the brightness parameter s 
as demonstrated in fig. 4. The relationship between the super-Poissonian character of the PCH and the molecular 
brightness e can be qualitatively understood. The average fluorescence intensity of a molecule in the excitation volume 
is characterized by the parameter 8. A particle with a larger value of e causes stronger intensity fluctuations as it enters 
and diffuses through the beam. The increase in the fluorescence intensity fluctuations leads to a further broadening 
of the PCH. This behavior is a consequence of the averaging of Poisson distributions over a wider intensity range as 
expressed by Mandel's formula. To quantify this statement, we define the fractional deviation Q, a measure of the 
deviation between the PCH and the Poisson distribution (Mandel, 1979), 




(32) 
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(Ak 2 )-{k) 



where 



and {k) are the variance and the expectation value of the photon counts, respectively, and y is the shape factor of the 
PSF (Thompson, 1991). The fractional deviation Q and the g{0) value of the intensity autocorrelation function (Thomp- 
son, 1991), 

(A/') (Ak')-(k) y 

s ^(iy--ir^=h (34) 

are closely related. The g(0) value is the ratio of the shape factor y to the average number of molecules N inside the 
PSF. The relationsh ip between the intensity moments and the factorial moments of the photon counts (Teich and Saleh , 
1988) allows us to express g(0) by the variance and the average of the photon counts. Thus, the ratio of Qto g(0) 
determines the average photon counts (k), according to eq. 1 0. 

[0099] A Poissonian distribution is defined by O=0, while super-Poissonian distributions require Q > 0 and sub- 
Poissonian distributions mandate Q < 0. The molecular brightness largely determines the super-Poissonian character 
of the PCH, since O is directly proportional to the parameter 8, which depends on the excitation power, the detection 
efficiency and the molecular species. 

Sensitivity of PCH algorithm 

[0100] To provide a quantitative description of the deviation of the photon counting histogram from a Poisson distri- 
bution, we define the following reduced Xp oj function, 



xlAe,K,M)= Mm 



z 

k=k tmn 



k -k 
*max ^rnin 



(35) 



where the standard deviation is given by. o k =(Mn(k,e,N)^-U(k,E,N))) y2 . The function % . describes the statistical 
significance of the deviation of the PCH from a Poisson distribution. If the value of % is fess than or equal to 1 then 

E-<?i 2 

the statistics are not sufficient to distinguish the data from a Poisson distribution. The larger the value of x poj! the 
stronger the deviation between the PCH and the Poisson distribution. Equation 35 requires the minimization of the 
function with respect to the parameter k . Since for most cases the value of k is nearly identical with the average photon 
counts of the PCH function, {k) = e/V, matching the first moments of the PCH and Poisson function is a good approx- 
imation of the function y 2 . 

[0101] We evaluated x poj for a variety of conditions to explore the sensitivity of PCH analysis. The dependence of 
Xp oj on the average number of molecules N for constant molecular brightness 8 is shown in fig. 5. Increasing the 
number of molecules at low concentrations results first in a steady increase of x pojl then the function reaches a max- 
imum and decreases monotonically at high particle concentrations. Thus, an optimal concentration exists, where the 
photon count distribution has a maximal deviation from the Poisson distribution. 

[0102] We can understand this result intuitively. Increasing the particle concentration leads to smaller amplitude 
fluctuations g(0), which means that the intensity distribution is approaching a delta function, as already mentioned 
before. Reducing the number of molecules in the observation volume produces stronger fluctuation amplitudes. How- 
ever, once the average number of molecules is less than one molecule, the probability that no molecule is found in the 
observation volume greatly increases. Therefore, except for the case where k= 0 the signal to noise ratio of the his- 
togram is markedly reduced. Thus, two effects shape the % . function and lead to a maximum around a particle con- 
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centration of one molecule in the observation volume. The function y . is shown in fig. 5 for two different molecular 

Poi 

brightness values to illustrate the influence of 8 upon the peak position. For brighter particles, the maximum of the 

function %p oj occurs at a lower concentration. An increase in the molecular brightness 8 increases the broadening of 

the PCH anci shifts the maximum of % 2 . to lower particle concentrations. 

5 [0103] The absolute values of the function depend strongly on the brightness parameter. In order to show the 

two curves on the same graph, the corresponding functions forthe larger brightness had to be scaled down to be visible 

(fig. 5). The dependence of yj . on the molecular brightness 8 was evaluated for two fixed concentrations (N = 0.01 

Poi 

and 1 .2) and both were fit to a straight line, thatyielded a decadic logarithm of about 1 .75 for both particle concentrations 
(fig- 6). 

10 [0104] The value of the x poj function is directly proportional to the number of data points M taken, and therefore 
directly proportional to the data acquisition time. If we formally define the square root of the Xp function to be a 
measure of the signal to noise ratio (SNR), then we find, that the SNR of PCH is proportional to the square root of the 
data acquisition time, and close to linearforthe molecular brightness. These conclusions are similarto those by Koppel, 
where the SNR was derived forthe autocorrelation function under the simplifying condition of Gaussian statistic (Koppel, 

15 1974). 

[0105] The dependence of the SNR on the number of molecules in the excitation volume is related to another study 
(Kask et al., 1 997). There, it was shown that the SNR of the second moment is constant at high particle concentrations, 
but decreases at low particle concentrations. However, the SNR of the third moment displays a maximum at a con- 
centration of about one particle per observation volume. This behavior of the SNR is also found in the SNR of the PCH, 
20 which is not surprising, since all moments are contained in the photon count distribution. 

PCH for multiple species 

[0106] Resolving multiple species is an important issue in many biological applications. Biological macromolecules 
25 interact with other molecules, and as a consequence of this network of interactions, the complex machinery of life is 
maintained. FCS has successfully been used to resolve multiple species (Rauer et al., 1 996). However, it is generally 
recognized that resolving two species by the autocorrelation function alone requires a difference in their diffusion co- 
efficients of the order of 2 or larger (Meseth et al., 1 999). This poses a severe restriction on the application of the pure 
autocorrelation approach to many biological systems, since the differences between the diffusion coefficients of bio- 
30 molecules is often less than a factor of two. The association of two monomer subunits, for example, to form a dimer 
is a widespread and important biological reaction mechanism (Berland et al., 1996). The increase in the diffusion 
coefficient for a dimer is about 25 percent of the monomer value. Since the diffusion coefficient approximately scales 
with the cubic root of the molecular mass, a difference in the molecular weight of about a factor of 8 would be required 
to resolve two components by the autocorrelation function alone. To address this intrinsic shortcoming of the autocor- 
35 relation approach two other methods have been introduced in the literature. One is based on higher order autocorre- 
lation functions (Palmer and Thompson, 1987; Palmer and Thompson, 1989) and another is based on higher order 
moment analysis (Qian and Elson, 1990a; Qian and Elson, 1990b). 

[0107] Here we introduce a new approach for separating fluorescent species based on the photon counting histo- 
gram. The histogram of photon counts is sensitive to the molecular brightness and was discussed in detail for a single 

40 species. If two species differ in their molecular brightness, then a molecule of the brighter species entering the obser- 
vation volume will produce a stronger fluorescent intensity change than the other species. By considering the statistics 
of these intensity changes, one can deduce the brightness and the average number of molecules of each species. 
Shot noise caused by the photodetection process is added to these intensity fluctuations. It was shown that the resulting 
photon counting statistics for multiple species is given by the consecutive convolution of the single species photon 

45 counting histograms. While the histogram of a single species requires two parameters, 2rparameters are required to 
describe the histogram for r species, namely the molecular brightness e,- and the average number of particles A/,for 
each species. In the following, we discuss the resolution of two species in detail, consider practical limitations and 
demonstrate the technique experimentally. 

50 Resolvability of two species 

[0108] We discuss the most difficult case where the two species must be resolved by the histogram alone without 
any further knowledge. From a practical point of view one wants to know what data acquisition time and concentrations 
to choose in order to resolve species of a given brightness. To address this question we calculated the theoretical 
55 histograms for different conditions, in order to identify experimentally favorable concentrations and brightness condi- 
tions. The theoretically determined two-species PCH function were then fit assuming a single species model and the 
reduced % 2 determined. A fit of a two species histogram by a single species model will result in a misfit, which gives 
rise to systematic residuals. The magnitude of the residuals tells us, whether it is feasible to distinguish between single 
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and multiple species. A reduced % 2 value of one or less indicates, that the data statistics is not sufficient to resolve the 
species, while a % 2 greater than one indicates that more than one species is present. We want to find out which bright- 
ness differences can be separated. The brightness for a given species is kept constant during the calculation of the 
histograms, but the concentration is varied in a systematic manner. For a fixed brightness ratio the results are best 
5 represented graphically in the form of a contour plot of the % 2 surface as a function of the logarithmic concentration of 
both species. The concentration of each species is expressed in number of molecules within the PSF. 
[0109] The above has been offered for illustrative purposes only, and is not intended to limit the scope of the invention 
of the application, which is as defined in the claims below. 
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Claims 

1 . A fluorescence spectrometer which comprises: 

a laser; at least one beam splitter positioned to receive a light beam from the laser and to divide said beam 
into a plurality of separate, first light beam portions; dichroic mirrors positioned to separately receive said first light 
beam portions and to reflect said beam portions at an angle to said first light beam portions; a plurality of transparent 
chambers for holding samples; objective lens systems respectively positioned in the path of said reflected beam 
portions to respectively focus each reflected beam portion to a point within one of said separate, transparent 
chambers; lenses respectively positioned to receive fluorescence from a sample for testing within said transparent 
chambers and to respectively focus said fluorescence at pinholes in opaque partitions, said lenses being positioned 
to receive said fluorescence which passes back through said objective lens system and said dichroic mirror; light 
detectors each respectively positioned adjacent to one of said partitions with one said partitions respectively po- 
sitioned between each of said lenses and said light detectors, to permit each light detector to sense fluorescence 
through said pinhole; and electronics to receive and process signals from each light detector. 

2. The spectrometer of Claim 1 in which a plate comprises said transparent chambers and a plurality of other trans- 
parent chambers for holding samples; and an x-y movement device which carries said plate and permits a first 
group of said chambers to be simultaneously respectively exposed to said first light beam portions for spectrometric 
analysis, followed by movement of the x-y device for exposure of another group of chambers to the first light beam 
portions for their spectrometric analysis. 

3. The spectrometer of Claim 2 in which a plurality of beam splitters divide said light beam into at least four first light 
beam portions which are spaced from each other for simultaneous irradiation of at least four of said transparent 
chambers, each separate first light beam portion irradiating only one of said transparent chambers. 

4. The spectrometer of Claim 3 in which said first light beam portions extend parallel to each other. 

5. The spectrometer of Claim 4 in which the chambers of said plate comprise open top wells, and the first light beam 
portions pass into said chambers through the bottoms of said wells. 

6. The spectrometer of Claim 5 in which fluorescence from each sample in each chamber passes in a straight path 
through said objective lens system, said dichroic mirror, said lens, and said pinhole to said light detector. 

7. The spectrometer of Claim 3 in which fluorescence from each sample in each chamber passes in a straight path 
through said objective lens system, said dichroic mirror, said lens, said pinhole to said light detector. 

8. The spectrometer of Claim 1 in which a plurality of dichroic mirrors divide said light beam into at least four first 
light beam portions which are spaced from each other for simultaneous irradiation of at least four transparent 
chambers, each separate first light beam portion irradiating only one of said transparent chambers. 

9. The spectrometer of Claim 8 in which said first light beam portions extend parallel to each other. 

10. The spectrometer of Claim 1 in which fluorescence from each sample in each chamber passes in a straight path 
through said objective lens system, said dichroic mirror, said lens, and said pinhole to said light detector. 

11. The method of performing fluorescence correlation spectroscopy, which comprises: 

acquiring a fluorescence measurement over time from a sample in the form of electronic raw data; electron- 
ically storing said raw data; and processing said raw data using a first algorithm without erasing said electronically 
stored raw data. 

1 2. The method of Claim 1 1 which comprises subsequently reprocessing said stored raw data using a second algorithm 
different from said first algorithm. 

13. The method of Claim 12 in which said raw data is acquired using a photon-mode technique. 

14. The method of Claim 11 in which said raw data is acquired using a photon-mode technique. 

15. The method of performing fluorescence correlation spectroscopy, which comprises: 
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acquiring a fluorescence measurement from a sample over time to obtain a histogram of photon counts; 
deriving the molecular brightness of at least one molecular species in said sample; and determining the concen- 
tration of said species. 

16. The method of Claim 15, in which the molecular brightness of a plurality of said species is determined from said 
fluorescence measurement. 

17. The method of Claim 15 in which the fluorescence measurement is taken from a sample volume of no more than 
1 0 femto liters. 

18. The method of Claim 17 in which the molecular brightness of a plurality of said species is determined from said 
fluorescence measurement. 

19. Fluorescence spectrometer which comprises: 

a laser; 

a dichroic mirror positioned to receive said light beam and to reflect a light beam portion at an angle to said 
light beam; 

a plate comprising a plurality of transparent chambers for holding samples, and an xy movement device which 
carries said plate and permits said chamber to be sequentially exposed to said light beam portion for spec- 
trometeric analysis, followed by movement of the x-y device for exposure of another chamber to the light beam 
portion for their spectrometeric analysis, an objective lens system positioned in the path of said beam portion 
to focus said beam portion to appoint within said separate transparent chamber; 

a lens system positioned to receive fluorescence from a sample for testing within said transparent chamber 
and to respectively focus said fluorescence at a pin hole in an opaque partition, said lens system being posi- 
tioned to receive said fluorescence which passes back through said objective lens system and said dichroic 
mirror; 

a light detector positioned adjacent to said partition, with said partition being positioned between said lens and 
said light detector, to permit said light detector to sense fluorescence through said pin hole; and 
electronics to receive and process signals from each light detector. 
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